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In order to directly compute physical two-pion K-decay amplitudes using lattice methods we must prepare a 
two-pion state with non-zero relative momentum. Building upon a proposal of Lellouch and Liischer, we describe 
a finite- volume method to realize such a state as the lowest energy state of two pions. 



1. Introduction 

The techniques of lattice gauge theory offer the 
possibility of accurate, non-perturbative calcula- 
tion of the masses and matrix elements predicted 
by low-energy QCD. Of particular interest are the 
properties of the K-meson system and the low- 
energy, weak matrix elements that are required to 
understand Kaon decays, especially those which 
enter the K — > tttt amplitudes which violate 
CP symmetry. Presently well-developed meth- 
ods permit the calculation of K to vacuum and 
K to 7T weak matrix elements. When combined 
with chiral perturbation theory and the accurate 
chiral symmetry of the domain wall fermion for- 
mulation, lattice methods can be used to evalu- 
ate all the matrix elements that contribute to the 
weak K — ► 7T7t amplitudes in the chiral limit |jl|,^) . 
In this paper we describe a method which may 
provide a practical approach to the direct lattice 
calculation of the desired, on-shell, K — ► ttw am- 
plitudes without the use of chiral perturbation 
theory. 

2. Finite volume decay amplitudes 

We begin by reviewing the finite-volume ap- 
proach developed by Lellouch and Liischer. To 
compute the matrix element of an operator O be- 
tween two physical states \A) and \B), one typi- 
cally evaluates the Euclidean Green's function: 

G(h,t 2 ,t 3 ) = Tr{B(t 1 )0(t 2 )A(t 3 )e- HT } 
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in the limit T 3> t\ 3> fa t 3 . Here the op- 
erators A and B are chosen to carry quantum 
numbers which insure that the desired states, \A) 
and \B) are each the lowest energy states that can 
be created from the vacuum by A and B, respec- 
tively. As demonstrated by Maiani and Testa||, 
this method fails when one of the desired states is 
a two-particle state above threshold. In this case 
the desired final state cannot be uniquely selected 
by its quantum numbers and the Euclidean time 
limit will give an unwanted, threshold amplitude. 

In a recent paper Q, Lellouch and Liischer pro- 
pose a method to circumvent this difficulty by 
tuning the finite lattice volume so that the first 
excited two-particle state allowed in that volume 
has the energy of the initial decaying particle. 
They work out an explicit formula which relates 
the resulting finite-volume matrix element to the 
desired infinite volume decay amplitude. The 
same result is obtained by different methods by 
Lin et ai, offering further insight^. As pro- 
posed, this method faces serious practical difficul- 
ties given our present computational capabilities. 
First, a large spatial volume of (6 Fermi) 3 is re- 
quired if the first-excited tt — tt state is to have 
an energy matching the Kaon mass. Second, the 
extraction of the next-leading exponential in the 
large-time limit is numerically difficult, especially 
if an accurate result is desired. Finally, while the 
matrix element of the 1 = 2 final state can be ob- 
tained from the coefficient of such a next-leading 
exponential, the more important 1 = state also 
receives a vacuum contribution implying that the 
desired amplitude is the third term in an expan- 
sion in increasingly small exponentials. 
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We propose to overcome these difficulties by 
imposing an anti-periodic boundary condition on 
the pions, say in the ^-direction, while using the 
usual periodic boundary conditions in the other 
two directions. With this choice the lowest en- 
ergy state of a single pion has momentum p z = |j 
instead of zero. With such a choice of boundary 
conditions, the lowest energy state with 1 = 2 
will be two pions, each with momentum p z w if;- 
Thus, the matrix element of such a state can be 
extracted from the leading exponential in the Eu- 
clidean Green's function. Since the vacuum state 
is unaffected by these boundary conditions, our 
1 = two-pion state will show the next-leading 
large-time behavior and a vacuum subtraction 
will be necessary for calculation of AI = i am- 
plitudes. 

The condition which insures that the two-pion 
state has the relative momentum required for a 
physical K — > tttt decay, now requires L z w 3 
Fermi, only one-half the lattice extent needed for 
the periodic case. The other two spatial dimen- 
sions are not constrained by these kinematics and 
might be chosen to be similar L x = L y = L z 
giving a reduced box size of (3 Fermi) 3 . Thus, 
the proposed method decreases the required spa- 
tial size from (6 Fermi) 3 to (3 Fermi) 3 and pro- 
motes the amplitude of interest to one of leading 
(/ = 2) or next-leading (/ = 0) large-time behav- 
ior. Since only periodic and cubically symmet- 
ric periodic boundary conditions were discussed 
in the original finite volume analysis || on which 
the Lellouch-Liischer paper Q is based, we must 
recompute Liischer's function <fi(k) for this new, 
less symmetrical case. This has been done with- 
out difficulty following the method presented in 
the original paper, Ref. Q. 

3. G-parity boundary condition 

It is not immediately obvious how to realize 
this boundary condition for the pion, because we 
have control only on the underlying quark fields. 
However, the G-parity operation gives the solu- 
tion since under G-parity: G\ir >= —\ir >. From 
this, we can see that anti-periodic boundary con- 
ditions for the pion can be achieved by performing 
a G-parity operation at the boundary — an oper- 



ation that also can be defined at the quark level: 
u — ► d and d — ► —u. (Recall that G is defined as 
the product of charge-conjugation and an isospin 
rotation, G = C expiirl y .) These unusual bound- 
ary were originally proposed and implemented by 
Wiese[Q and similar boundary conditions have 
also been used in Ref. ||. 

Since we must also treat the Kaon, we must 
impose a charge conjugate boundary condition on 
the strange quark for consistency. This is easily 
done by including a fictitious charmed quark, de- 
generate with the strange quark and extending 
the usual G-parity by treating the (c, s) quark 
pair as an independent iso-doublet. Our initial 
finite volume eigenstate can then be constructed 
as a G-even mixture (K° + D°) j\pl. However, be- 
cause our weak decay operator and our final states 
have definite strangeness and charm, this mixing 
only introduces an extra "finite-volume factor" of 
\j\[2 which can easily be removed. 

4. Lattice Action 

4.1. Symmetry 

By imposing this G-parity boundary condition, 
we break chiral symmetry in a manner similar 
to the chiral symmetry breaking of a mass term. 
This was studied for the continuum in Ref. Q|. 
On the lattice, the boundary condition is em- 
bedded in the action itself, adding a boundary- 
specific term. For example, with G-parity bound- 
ary conditions we have a term in the action den- 
sity of the form: 

u(L z - l)U 3 (L z - l) 73 Cd(0) T , (1) 

for Dirac conventions in which charge conjugation 
is represented by q —> Cq T . Here u(L z — 1) is eval- 
uated at the site z = L z — 1 and d(0) at the site 
z = while U 3 (L Z — 1) is the link variable join- 
ing these two sites across the boundary. Under 
chiral SU(2)l ® SU(2)r such a term transforms 
as a (5,5). While such a boundary term may 
appear to break translational invariance, the fact 
that the lattice Lagrangian is invariant under the 
G-parity operation ensures us that the system has 
translational symmetry, if appropriately defined. 
(Note, when exploring unconventional boundary 
conditions it is important to preserve the conser- 
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vation of lattice momentum and hence to require 
translational invariance.) 

4.2. Implementation 

The actual implementation of these boundary 
conditions represents a fairly minor modification 
of the usual lattice Dirac operator. An easy way 
of visualizing these G-parity boundary conditions 
begins with doubling the lattice in the z-direction, 
with the complex conjugates of gauge links in the 
original volume copied onto the links translated 
by L z into this extended volume. One then eas- 
ily defines a single-flavor Dirac operator on this 
L x x L y x 2 • L z lattice in the usual way, but 
imposing anti-periodic boundary conditions be- 
tween the z = and z = 2L Z — 1 boundaries. If 
a quark field in the original volume is identified 
as a u quark, then the quark field in the doubled, 
L z < z < 2L Z — 1 volume would be viewed as a d 
quark. Note, the computational load is twice that 
of a usual calculation on a L x x L y x L z volume — 
our boundary conditions have effectively required 
separate treatment for the u and d quarks. Such 
an approach works for both Wilson and domain 
wall fermions. 



5. Conclusion 

Since the G-parity boundary condition is im- 
posed in only one direction, we lose cubic sym- 
metry. In the cubically symmetric case the state 
of interest with orbital angular momenta I = is 
mixed only with states having I differing by mul- 
tiples of 4. However, without this cubic symme- 
try only parity invariance remains which permits 
mixing with all even I. Thus, in the practical 
application of the Lellouch and Liischer method 
where the mixed, higher-/ states are treated as 
free and their phase shifts ignored, we must ne- 
glect the I = 2 phase shift 8^ (mj{), This reduces 
the accuracy of this approximation, introducing 
errors of order (R/L) 4 rather than (R/L) s as is 
the case for cubic symmetry Here R is a typi- 
cal strong interaction distance, e.g. R = l/m p . 

There is an additional complexity for the 1 = 
final state that should be discussed. For simplic- 
ity, we consider a full QCD calculation which in- 
cludes only dynamical u and d quarks. For such a 



system, in addition to the three pions, there will 
be a fourth SU{2) singlet, quark-anti-quark state 
called here the rf , This rj state is G-parity even 
in contrast to G-parity odd triplet of pions and 
hence will obey even boundary conditions. Thus 
an intended two-pion, 1 = state may mix with 
an unwanted, threshold two-?/ state which will 
have the same quantum numbers. Fortunately, 
in the two-flavor calculation just described, ax- 
ial anomaly effects are expected make the non- 
Goldstone rj heavier than the physical Kaon so 
this two-?/ state will play no role. However, in 
a quenched calculation, this two-?/ state will be 
essentially degenerate with the desired 1 = two- 
pion state and obstruct the calculation. In fact, 
the hairpin diagrams entering the quenched two- 
rj' amplitudes will give an extra power of t or t 2 
further obscuring the state of interest. Thus, this 
G-parity boundary condition can be used in the 
quenched approximation only for the evaluation 
of Al = 3/2 amplitudes. 

In summary, we propose imposing anti-periodic 
boundary conditions on the pion by implement- 
ing a G-parity operation at boundary. This 
reduces the required simulation volume to w 
(3 Fermi) 3 and, more importantly, allows the 
physical {^^(I^K^^K) amplitude to be ex- 
tracted from the leading/next-leading large-time 
behavior of a Euclidean correlation function for 
the 1 = 2/1 = final states. We hope to try a 
calculation according to this method in the com- 
ing year and thank our colleagues in the RBC col- 
laboration for useful suggestions and criticisms. 

REFERENCES 



J. I. Noaki, et a l |hcp-lat/0108013 



T. Blum, et al, |hep-lat/0110075 . 
L. Maiani and M. Testa, Phys. Lett. B245 
(1990) 585-590. 

L. Lellouch and M. Luscher, Comm. Math. 
Phys. 219 (2001) 31-44. 

C. J. D. Lin, et al, Nucl.Phys. B619 (2001) 
467-498. 

M. Luscher, Nucl. Phys. B354 (1991) 531-578. 
U. J. Wiese, Nucl. Phys. B375 (1992) 45-66. 
J. M. Carmona, et al, Int. J. Mod. Phys. C 
11, 637 (2000). 



